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ABSTRACT 

We present strong observational evidence for a relationship between the direction of 
a pulsar's motion and its rotation axis. We show carefully calibrated polarization 
data for 25 pulsars, 20 of which display linearly polarized emission from the pulse 
longitude at closest approach to the magnetic pole. Such data allow determination of 
the position angle of the linear polarisation which in turn reflects the position angle 
of the rotation axis. Of these 20 pulsars, 10 show an offset between the velocity vector 
and the polarisation position angle which is either less than 10° or more than 80°, 
a fraction which is very unlikely by random chance. We believe that the bimodal 
nature of the distribution arises from the presence of orthogonal polarisation modes 
in the pulsar radio emission. In some cases this orthogonal ambiguity is resolved by 
observations at other wavelengths so that we conclude that the velocity vector and 
the rotation axis are aligned at birth. Strengthening the case is the fact that 4 of the 5 
pulsars with ages less than 3 Myr show this relationship, including the Vela pulsar. We 
discuss the implications of these findings in the context of the Spruit & Phinney (1998) 
model of pulsar birth-kicks. We point out that, contrary to claims in the literature, 
observations of double neutron star systems do not rule out aligned kick models and 
describe a possible observational test involving the double pulsar system. 
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1 INTRODUCTION 

The velocities of pulsars are significantly larger than those 
of their progenitor (high-mass) stars (e.g. Lyne & Lorimer 
1994). This implies that the birth process of pulsars also 
produces their high velocities and that the supernova or 
events soon thereafter must be asymmetrical. The mecha- 
nisms driving this asymmetry are far from clear. Tademaru 
& Harrison (1975) proposed the so-called rocket mechanism 
whereby an offset of the magnetic dipole from the centre 
of the star causes the star to accelerate along its rotation 
axis. Observational tests of this mechanism by Morris, Rad- 
hakrishnan & Shukre (1976) and Tademaru (1977) reached 
opposing conclusions, with the former claiming no correla- 
tion between the velocity vector and the spin axis, and the 
latter claiming a good correlation. Following this, Anderson 
& Lyne (1983) revisited the issue with an updated database 
of pulsar proper motions. They also saw no apparent align- 
ment between the proper motion angle and the rotation axis 
although they included many old pulsars in their sample. 



whose direction of motion may have significantly evolved 
in the Galactic potential. In the late 1980s, a group of pa- 
pers concentrated mainly on asymmetric explosions and the 
break-up of binary systems as the origin for pulsar veloci- 
ties (Radhakrishnan & Shukre 1986; Dewey & Cordes 1987; 
Bailes 1989) and the rocket model fell out of favour. 

A decade later, the debate started anew when Spruit 
& Phinney (1998) and Cowsik (1998) proposed a model in 
which the same (off-centre) explosion processes that give 
pulsars their high velocity also produce their fast rotation. 
Although their models were motivated by theoretical results 
suggesting that the pre-supernova stellar configuration is 
not spinning quickly enough to account for birth periods 
of ~10 ms from conservation of angular momentum alone, 
the basic physics applies - any significant off-centre impulses 
would both accelerate and torque the pulsar. In the Spruit 
& Phinney (1998) model, a single, short-duration impulse 
leads to the rotation axis being perpendicular to the veloc- 
ity vector. However, even a few short duration kicks as the 
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star spins serve to destroy this correlation. They also consid- 
ered kicks of finite duration (> the spin period) and showed 
in this case that the rotation axis should be parallel with 
the velocity vector. Deshpande, Ramachandran & Radhakr- 
ishnan (f999) tested these ideas by using the literature to 
search for pulsars which had accurate proper motion mea- 
surements and good polarization observations. They used 
the polarization measurements to determine the point of in- 
flexion of the PA sweep although give no details as to how 
this was done. They took the associated position angle as 
the angle of the rotation axis on the sky (or, to allow for 
orthogonal mode emission, its orthogonal). In their sample 
of 29 pulsars they saw no statistical evidence for any cor- 
relation between the velocity vector and the rotation axis. 
In the context of the Spruit & Phinney (1998) model they 
favoured multiple short duration impulses as the origin of 
the spin and velocity kicks. 

Observational evidence for a correlation between the 
spin axis and the proper motion next came from an unex- 
pected source. Helfand et al. (2001) showed high resolution 
X-ray images of the Vela pulsar and its surrounding nebula. 
Two prominent arcs are seen, which Helfand et al. (2001) 
interpret as originating from the front surface of an equato- 
rial wind whose axis of symmetry is the rotation axis of the 
pulsar. The position angle of the axis of this torus is 130° 
and 310° . The axis of the torus, created by pulsar winds, is 
likely aligned with the pulsar rotation axis. Helfand et al. 
(2001) then note that the position angles of the torus and 
pulsar spin axes are within 10° of the proper motion vector 
of 301° obtained by Dodson et al. (2003). Following from 
the Vela result, X-ray imaging of a number of young pulsars 
showed that many had pulsar wind nebulae consisting of an 
equatorial torus and polar jets. Ng & Romani (2004) used 
this information to determine the angle of rotation axis on 
the plane of they sky and compared this to the velocity vec- 
tor if available. For the Crab and Vela pulsars these clearly 
matched but there was also good evidence for a match in the 
other pulsars in their sample (PSRs J0538-I-2817, B1951+32 
and B 1706-44). It appears likely from this study that the 
directions of motion of young pulsars lie along the rotation 
axes, although their sample is small. 

Historically, the velocity vectors of pulsars have been 
measured either through proper motion measurements us- 
ing Very Long Baseline Interferometry (VLBI) techniques 
or through highly accurate timing. The first method tends 
to be possible only for relatively nearby pulsars, whereas the 
second method weis generally restricted to millisecond pul- 
sars which show almost no intrinsic timing jitter. Recently, 
however, Hobbs et al. (2004) have shown that it is possi- 
ble to measure proper motions of pulsars even in the pres- 
ence of timing noise. In brief, this involves very long (tens 
of years) timing sets coupled with a 'whitening' technique 
to remove the timing noise. Using this method, Hobbs et 
al. (2004) list proper motion measurements for 302 pulsars. 
Updates to these parameters and a catalogue of the best 
available proper motions from cither timing or interferom- 
etry are available in Hobbs et al. (2005) and in particular 
provide celestial position angles (PA„) of the velocity vec- 
tors. 

In the absence of external features such as pulsar wind 
nebulae and tori, the position angle of the axis of rotation 
(PA,.) of a pulsar has to be determined from polarisation 




Figure 1. Geometry of the pulsar beam and velocity vector as 
seen by tlie observer. Tlie observer sees a projection of tlie di- 
rection of the rotation axis, PA^, and the velocity vector, PA„, 
onto the plane-of-the sky shown on the right. The angle a is the 
inclination angle of the magnetic and rotation axis, while /3 is 
the impact parameter between the magnetic axis and the line-of- 
sight at closest approach. The vector PAy lies along the field line 
in this fiducial plane ((/>o) and hence has the same value as PA^. 
PAx lies orthogonal to the magnetic field at closest approach. 
In the absence of other information we do not know whether the 
observed PAo corresponds to PAy or PAx- Note that PAy and 
PAx (and hence PAr) have 180° ambiguities inherent to all po- 
larization measurements, unlike PA^. 

observations using a several-step procedure. One first has 
to determine PAo, the position angle of polarization at the 
point of closest approach of the observer's line of sight to the 
magnetic pole. This angle must then be corrected for Fara- 
day rotation in the interstellar medium and the ionosphere 
in order to reflect the true angle at the pulsar. Determina- 
tion of PAo itself is difficult. We discuss two techniques. 

In the rotating vector model (RVM) of Radhakrishnan 
& Cooke (1969), the radiation is beamed along the field lines 
and the plane of polarization is determined by the angle of 
the magnetic field as it sweeps past the line of sight. The 
PA as a function of pulse longitude, (f), can be expressed as 

/ sinasin(c6 — 0o) \/,^ 

PA = PAo + axctan — ^ \ . . . 1 

y smQ cosa — cosC sma cos(^ — q>o) J 

Here, a is the angle between the rotation axis and the mag- 
netic axis, where C, = a + fi being the angle at closest 
approach of the line of sight to the magnetic axis. 00 is the 
corresponding pulse longitude at which the PA is then PAq. 
Unfortunately, for most pulsars it is difficult to determine a 
and (3 with any degree of accuracy, partly because the longi- 
tude over which pulsars emit is rather small and partly be- 
cause strong deviations from a simple swing of PA are often 
observed. This makes the determination of PAo straightfor- 
ward only in the ~15 per cent of pulsars for which the RVM 
works (see also Everett & Weisberg 2001). 

It is possible, however, to determine 4>o and hence PAo 
without RVM fitting. The frequency evolution of pulsar pro- 
files has shown that components close to (^o have steeper 
spectra than components at larger values of (Rankin 1983; 
Lyne & Manchester 1988; Kramer 1994). Also, there is of- 
ten change in the handedness of circular polarization at 0o 
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Table 1. Proper motions, rotation measures and polarisation position angles for a sample of 25 pulsars. 



Age Dist Proper Motion Rotation Measure Poln 

.liiaiiie Biiaiiic log(r,.) fia Ref PA,. Previous Ref Measured PAo 

(yr) (kpc) l_ma^/\-r) (inas/vr) (iI<'k) Israel ui ") (rad m ") (deK) 



J0525+1115 


B0523+11 


7.9 


3.1 


27(2) 


-24(10) 


1 


132(16) 


35(3) 


1 


37(2) 


-65(4) 


J0601-0527 


B0559-05 


6.7 


3.9 


-5(3) 


-20(6) 


1 


194(16) 


67(5) 


2 


62.4(20) 


— 


J0630-2834 


B0628-28 


6.4 


1.4 


-44.6(9) 


19.5(22) 


2 


294(3) 


46.19(9) 


3 


46.53(12) 


26(2) 


J0742-2822 


B0740-28 


5.2 


2.1 


-29(2) 


4(2) 


3 


278(5) 


156(5) 


4 


149.95(5) 


-81.7(1) 


J0820-1350 


B0818-13 


7.0 


2.0 


18(3) 


-46(5) 


1 


159(6) 


-1.2(4) 


3 


-7.2(12) 


65(2) 


J0835-4510 


B0833-45 


4.1 


0.29 


-49.60(6) 


29.8(1) 


4 


301.0(1) 


38.17(6) 


5 


31.38(1) 


36.8(1) 


J0908-1739 


B0906-17 


7.0 


0.9 


22(3) 


-94(4) 


1 


167(2) 


-36(6) 


2 


-31(4) 


— 


J0922+0638 


B0919+06 


5.7 


1.2 


18.35(6) 


86.56(12) 


5 


12.0(1) 


32(2) 


2 


29.2(3) 


— 


J0953+0755 


B0950+08 


7.2 


0.26 


-2.09(8) 


29.46(7) 


6 


355.9(2) 


1.35(15) 


3 


-0.66(4) 


14.9(1) 


J1136+1551 


B1133+16 


6.7 


0.36 


-74.0(4) 


368.1(3) 


6 


348.6(1) 


3.9(2) 


6 


1.1(2) 


-78(2) 


J1239+2453 


B1237+25 


7.4 


0.86 


-106.82(17) 


49.92(18) 


6 


295.0(1) 


-0.33(6) 


3 


-2.6(4) 


-66(1) 


J1430-6623 


B1426-66 


6.7 


1.0 


-31(5) 


-21(3) 


7 


236(9) 


-12(3) 


3 


-19.2(3) 


-28.5(7) 


J1453-6413 


B1449-64 


6.0 


2.1 


-16(1) 


-21.3(8) 


7 


217(3) 


-22.3(14) 


3 


-18.6(2) 


-56.9(4) 


J1456-6843 


B1451-68 


7.6 


0.45 


-39.5(4) 


-12.3(3) 


8 


252.7(6) 


-2(3) 


3 


-4.0(3) 


-31.6(6) 


J1645-0317 


B1642-03 


6.5 


1.1 


-3.7(15) 


30.0(16) 


2 


353(3) 


15.8(3) 


3 


17(2) 


56(4) 


J1709-1640 


B1706-16 


6.2 


0.8 


-9(3) 


-44(21) 


1 


192(16) 


-1.3(3) 


2 


-2(5) 


15(2) 


J1709-4429 


B1706-44 


4.2 


2.3 








160(10) 


-7(4) 


7 


0.70(7) 




J1740+1311 


B1737+13 


6.9 


1.5 


-21.5(22) 


-19.7(22) 


2 


228(6) 


64.4(16) 


1 


65(2) 


-46(4) 


J1844+1454 


B1842+14 


6.5 


2.2 


18(4) 


25(7) 


1 


36(15) 


121(8) 


2 


109.0(13) 


-52(2) 


J1900-2600 


B1857-26 


7.7 


2.0 


-19.9(3) 


-47.3(9) 


9 


202.8(7) 


-7.3(8) 


3 


-2.3(8) 


-43(2) 


J1913-0440 


B1911-04 


6.5 


2.8 


6(2) 


-24(7) 


1 


166(11) 


12(3) 


2 


4.4(9) 


-68(2) 


J1921+2153 


B1919+21 


7.2 


1.1 


19(4) 


28(6) 


1 


34(12) 


-16.5(5) 


3 


-13(3) 


-35(2) 


J1932+1059 


B1929+10 


6.5 


0.36 


94.03(14) 


43.4(3) 


5 


65.2(2) 


-6.1(1.0) 


3 


-6.87(2) 


-11.3(1) 


J1935+1616 


B1933+16 


6.0 


5.6 


1.1(4) 


-17.6(7) 


1 


176.4(14) 


-1.9(4) 


6 


-10.2(3) 


10.1(7) 


J1941-2602 


B1937-26 


6.8 


1.7 


12(2) 


-10(4) 


2 


130(17) 


-41(4) 


7 


-33.5(8) 





Proper motion references: 1. updated parameters since Hobbs et al. (2004), 2. Brisken et al. (2003), 3. Fomalont et al. (1997), 4. 
Dodson et al. (2003), 5. Chatterjee et al. (2001), 6. Brisken et al. (2002), 7. Bailes et al. (1990b), 8. Bailes et al. (1990a), 9. Fomalont et 
al. (1999). Rotation measure references: 1. Weisberg et al. (2004), 2. Hamilton & Lyne (1987), 3. Hamilton, McCulloch & Manchester 
(1981; unpublished), 4. Han, Manchester & Qiao (1999), 5. Hamilton et al. (1987), 6. Manchester (1972), 7. Qiao et al. (1995) 



(Radhakrishnan & Rankin 1990). Finally, symmetry in pulse 
profiles is also an important indicator of (f>o (Rankin 1983; 
Lyne & Manchester 1988). Most of the pulsars in our sample 
have been extensively observed at a wide range of frequen- 
cies. We use as much information as possible to locate (bo 
in the polarization profiles with the most important being 
RVM fits where robust. 

An added complication is that pulsar emission can occur 
in two orthogonal modes which have PAs at right angles to 
each other (Manchester, Taylor & Huguenin 1975; Backer, 
Rankin & Campbell 1975). This can distort the swing of 
PA with pulse longitude but also provides a 90° ambiguity 
in any analysis of polarized alignment as already noted by 
Tademaru (1977). Figure 1 illustrates the angles in the prob- 
lem and shows the orthogonal emission possibilities. It can 
be seen that although one would like to associate the mea- 
sured PAo directly with PA,., it may be that PA^ = PAo + 
90° if the pulsar emission is in the orthogonal mode. 



2 MOTIVATION AND SOURCE SELECTION 

Given the recent results, the time seems right to re- visit the 
question, from an observational point of view, of whether 
pulsar spin axes and velocity vectors align. In order to 
achieve this, we need proper motion measurements with 



small errors, high-quality polarization data with absolute 
position angles, and accurate rotation measures (RM) to re- 
move Faraday rotation from the measured polarization PAs. 
Our intention was to make high-quality polarization obser- 
vations on a number of pulsars for which good proper motion 
data exist. 

There were three main criteria for source selection. 
First, the chosen pulsars had to have accurate proper mo- 
tions, with an angle on the sky known to better than 20°. 
Secondly, the pulsars had to be relatively young so that the 
velocity vector is an accurate representation of their birth 
velocity vector (i.e. had not had time to bo significantly con- 
taminated by acceleration in the Galactic potential) . Finally, 
we concentrate here on southern pulsars with declinations 
less than -|-25°; a future study of northern pulsars is planned. 
We chose a total of 25 pulsars which met most of these crite- 
ria. We have included PSR J1709— 4429 in our sample oven 
though it does not have a measured proper motion. Dod- 
son & Golap (2002) and Bock & Gvaramadze (2002) have 
proposed that the pulsar is associated with the supernova 
remnant G343. 1—2.3 and claim that the proper motion vec- 
tor should thus be in the range 150° to 170°. Although the 
claimed association has been disputed (Nicastro, Johnston 
& Koribalski 1996; Frail & Scharringhausen 1997), Ng & 
Romani (2004) derived a position angle vector of 175° for 
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the rotation axis based on observations of an X-ray torus 
around the pulsar. 

Table 1 lists the pulsars, none of which are binary or 
millisecond pulsars. The first two columns list the J and B 
names of the pulsars. Column 3 gives the logarithm of the 
pulsar's characteristic age, Tc = P/2P, with P the pulsar 
period and P the period derivative. Column 4 gives the dis- 
tance to the pulsar in kpc. Columns 5 to 8 give the proper 
motions in RA and Dec, the reference and PA„, the position 
angle of the velocity vector on the sky measured counter- 
clockwise from north. In all cases the number in brackets 
indicates the error on the last digit (s). There are 10 pulsars 
in common with the Dcshpandc ct al. (1999) sample and 5 
in common with Anderson & Lyrie (1983). In most cases, the 
errors in the direction of the velocity are now significantly 
reduced. 

In the following sections we describe the polarization 

observations and the procedures required to determine ab- 
solute position angles at the pulsar, including the necessary 
accurate RM determination. We present our results in sec- 
tion 5 and 6 and discuss the implications in section 7. 



3 OBSERVATIONS 

The observations were caxried out using the Parkes radio 
telescope from 2004 November 27 to 30. For the first two 

days we used the H-OH receiver at a central frequency of 
1369 MHz with a bandwidth of 256 MHz. On the last two 
days we used the 10/50 cm receiver at a central frequency 
of 3100 MHz with a bandwidth of 512 MHz. Both receivers 
have orthogonal linear feeds and also have a pulsed calibra- 
tion signal which can be injected at a position angle of 45° 
to the two feed probes. A digital correlator was used which 
subdivided the bandwidth into 1024 frequency channels and 
provided all four Stokes' parameters. We also recorded 1024 
phase bins per pulse period for each Stokes' parameter. 

The pulsars were observed twice for 20 mirmtes each 
time, with a feed rotation of 90° between observations. Prior 
to the observation of the pulsar a 3-min observation of the 
pulsed calibration signal was made. The data were writ- 
ten to disk in FITS format for subsequent off-line analysis. 
Data analysis was carried out using the PSRCHIVE software 
package (Hotan, van Straten & Manchester 2004). 



4 DATA ANALYSIS AND CALIBRATION 

In order to determine the intrinsic positional angle of lin- 
ear polarization at the pulsar, there are a number of fac- 
tors which need to be taken into account. Working back- 
wards from the backend system to the pulsar, these are (i) 
the phase delay in the paths from the feed to the backend, 
(ii) gain differences between the two orthogonal probes of 
the feed, (iii) possible cross-talk or leakage between the two 
probes, (iv) the orientation of the feed probes with respect 
to the meridian, (v) the parallactic angle of the source, (vi) 
the rotation measure due to the terrestrial ionosphere and 
(vii) the rotation measure due to the interstellar medium. 

The first two of these can be removed using standard 
calibration techniques. We observed a pulsed calibration sig- 
nal, injected at 45° to the two probes so that the signal is 



100 per cent polarised in Stokes U. Gain corrections can be 
made by comparing Stokes Q with Stokes / and phase cor- 
rections can be made by comparing Stokes U with Stokes V 
on a channel by channel basis in the output data. Cross talk, 
or leakage terms can be largely eliminated by summing the 
two scans made with 90° feed offsets (Johnston 2002). As 
a further check we observed PSR J1359-6038 on 8 separate 
occasions over a range of hour angles. This allowed us to 
determine impurities in the feed, which we estimate to be 
no greater than 2 percent. 

The RM must be determined to high accuracy because 
of the need to account for the Faraday rotation of the po- 
larized signal. If the error in the RM is denoted S-rm then 
the error Spa (in radians) in the Faraday-corrected PA, can 
be calculated as 



where c is the speed of light and v the observing frequency. 
At 1369 MHz, therefore, an RM error of 1 rad m~^ yields 
a PA error of 2.8°. It has been noted by a number of au- 
thors (e.g. Weisberg et al. 2004) that RMs vary with epoch 
and we thus measured RMs for all the pulsars in our sam- 
ple rather than use previous values which can be decades 
old. There arc two ways to determine the RM. The first is 
to measure the PA rotation across the 256 MHz bandwidth 
at 1369 MHz. The second method is to determine the dif- 
ferences in PA obtained at 1369 and 3100 MHz. The latter 
gives more accurate results in principle but profile shapes 
and the intrinsic PA can sometimes be frequency-dependent 
(Karastergiou & Johnston 2005). For all the pulsars in our 
sample for which we have dual-frequency observations, we 
measured the RM using the latter technique. For pulsars 
which we only observed at 1369 MHz, we measured the RM 
across the band. 

The contributions to the measured values of RM arising 
in the ionosphere have been estimated by integrating a time- 
dependent model of the ionospheric electron density and ge- 
omagnetic field through the ionosphere along the sight line 
between the telescope and the pulsar using code provided 
by JL Han. Values of the ionospheric RM range from to 
—2 rad m~^. We have subtracted the ionospheric RMs from 
the measured values to leave only the interstellar compo- 
nents of the RM. 

As a final chock of the entire calibration technique, we 
also carried out short observations of PSR J 1359-6038 at 
1384 MHz with the Australia Telescope Compact Array 
(ATCA), a 6 element interferometer with a correlator which 
is capable of (crude) pulsar binning. We chose PSR J1359- 
6038 because it has a flat PA swing across the pulse and 
is not therefore affected by time averaging. We determined 
the absolute PA of the radiation using the standard package 
MIRIAD. The PA at the ATCA was the same as that mea- 
sured at Parkes within the errors, after the RM was taken 
into account. This gives us confidence that our calibration 
techniques at Parkes are working effectively and that we 
have correctly obtained absolute PAs. 

5 RESULTS 

The measured values of interstellar RM are listed in columns 
9 to 11 of table 1 as well as previous measurements, along 
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Figure 2. Polarization profile of the Vela pulsar at 1.4 GHz. The 
bottom panel shows the total intensity (black, thick line), linear 
polarization (grey line) and circular polarization (pale grey, thin 
line). The top panel shows the position angle of the polarized 
radiation as a function of pulse longitude. Two curves are shown. 
The first (dark points) is the measured PA, the second (lighter 
points) is the measured PA + 90°. The position angle of the 
velocity vector, PA„, is shown as a solid, horizontal line. Pulse 
longitude 0.0 ((/>o) is the closest approach of the line of sight to 
the magnetic pole as determined from the rotating vector model 
fit. The PA at this longitude is PAq from Equation 1. See also 
Fig. 1 for a depiction of the relevant angles. 



with the corresponding reference. The error in RM is typi- 
cally less than 0.5 rad for pulsars with dual frequency 
measurements but about 10 times larger for the other pul- 
sars. The small error in RM is crucial for determining abso- 
lute position angles. 

There are a number of cases in which the RM differs 
significantly from previous values, in particular that of the 
Vela pulsar. In 1977, the RM of the Vela pulsar was found 
to be 38 rad m~^ (Hamilton et al. 1977b). In a later paper, 
Hamilton et al. (1985) showed that the RM had increased 
to 46 rad m~^ and in fact appeared to be linearly increas- 
ing since at least 1970. As far as we are aware there are 
no further RM measurements quoted in the literature. Our 
measured value of the RM (2004 November) shows it to have 
a value of 31.4 rad m~^, significantly different to the older 
measurements. Clearly, following an RM increase for at least 
15 years (from 1970 to 1985), the RM has now significantly 
decreased and is currently below the 1970 value. At the time 
where the RM was increasing, the dispersion measure (DM) 
was decreasing (see Hamilton et al. 1985). We are not aware 
of any further DM measurements for Vela but note that DM 



values depend strongly on profile alignment between differ- 
ent frequencies and care should be taken when comparing 
results from different epochs. 

In the rest of this section, we present polarisation mea- 
surements of the individual sources and describe the deter- 
mination of PAo, the intrinsic position angle of polarization 
at 4>o, the longitude of closest approach of the line-of-sight 
to the magnetic pole. Clearly, our choice of 00 is crucial to 
our arguments and we justify this choice in detail for each of 
the pulsars in the paragraphs below. In particular we have 
taken into account evidence such as profile symmetry, profile 
frequency evolution and the circular polarization in deter- 
mining our choice, in a similar fashion to the morphology 
work of Rankin (1983) and Lyne & Manchester (1988). The 
final column of table 1 contains the measured value of PAq 
and its error which is derived from the error in the RM de- 
termination for each pulsar. Apart from the choice of 0o 
itself, additional observational errors due to signal to noise 
issues contribute only a very small amount. We deal with 
other sources of error in Section 6.3 below. 

Figures 2 and 3 show the polarization profiles for the 
25 pulsars in the sample. We first present the results for the 
Vela pulsar as an illustrative case study, and then the results 
for the other pulsars in order of increasing Right Ascension. 
PSR J0835-4510 (Vela; Figure 2): Figure 2 shows the 
polarization profile of the Vela pulsar from our data. This 
pulsar represents one of the rare cases where an RVM fit is 
robust, so we choose 00 from the fit. The RVM fit constrains 
both a and /3 reasonably well (—137° and 6.5° for these data 
using the sign conventions in Everett & Weisberg 2001). The 
fitted location of the closest approach of the line of sight to 
the magnetic pole is at longitude 0.0 in the figure. 
PSR J0525+1115 (Figure 3a): The profile at 1.4 GHz 
consists of at least three components. The linear polarization 
is low throughout and the circular polarization changes sign 
between the first and second component. Following Weisberg 
et al. (1999), we identify 0o to be at the longitude of the 
circular polarization sign change. 

PSR J0601-0527 (Figure 3b): At this frequency the PA 
consists of two flat sections, near —42° in the leading compo- 
nent and at -f 75° in the trailing component. We have placed 
00 where the sign of circular polarization changes. At this 
location there is no linear polarization and hence we cannot 
measure PAo. It is possible that 0o is actually somewhat 
earlier as the leading component is much more prominent in 
low frequency observations (Gould & Lyne 1998). Further- 
more, at low frequencies there does appear to be a swing of 
PA which connects the two fiat regions although this is com- 
plicated by an orthogonal jump. Absolute PA determination 
at a lower frequency may resolve this issue. 
PSR J0630-2834 (Figure 3c): The profile is rather sym- 
metric in this pulsar and we have placed 0o at the centre of 
symmetry. There is also rather little evolution of the pulse 
shape with frequency and it seems likely that the entire emis- 
sion originates from a component near 0o. 
PSR J0742-2822 (Figure 3d): As in the case of the Vela 
pulsar, the PA swing of this pulsar permits a robust RVM fit 
to be obtained and we used the fitted value of 0o to locate 
the closest approach to the magnetic pole. This longitude is 
also the location of the mid-point of the profile at the 10 per 
cent emission level. 

PSR J0820-1350 (Figure 3e): The pulsar shows rather 
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Figure 3. (a)-(d). Polarization profiles at 1.4 GHz for PSR J0525+1115 (top left), PSR J0601-0527 (top right), PSR J0630-2834 
(bottom left) and PSR J0742— 2822 (bottom right). See caption for Fig. 2 and text for details. The choice of longitude zero for each 
pulsar is given in the text and is marked with a dashed line where an accurate determination can be made. 
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Figure 3. (e)-(h). Polarization profiles at 1.4 GHz for PSR J0820-1350 (top left), PSR J0908-1739 (top right), PSR J0922+0638 
(bottom left) and PSR J0953-I-0755 (bottom right). See caption for Fig. 2 and text for details. The choice of longitude zero for each 
pulsar is given in the text and is marked with a dashed line where an accurate determination can be made. 
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Figure 3. (i)-(l). Polarization profiles at 1.4 GHz for PSR J1136+1551 (top left), PSR J1239+2453 (top right), PSR J1430-6623 (bottom 
left) and PSR J1453— 6413 (bottom right). See caption for Fig. 2 and text for details. The choice of longitude zero for each pulsar is 
given in the text and is marked with a dashed line where an accurate determination can be made. 
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Figure 3. (m)-(p). Polarization profiles at 1.4 GHz for PSR J1456-6843 (top left), PSR J1645-0317 (top right), PSR J1709-1640 
(bottom left) and PSR J1709— 4429 (bottom right). See caption for Fig. 2 and text for details. The choice of longitude zero for each 
pulsar is given in the text and is marked with a dashed line where an accurate determination can be made. 
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Figure 3. (q)-(t). Polarization profiles at 1.4 GHz for PSR J1740+1311 (top left), PSR J1844+1454 (top right), PSR J1900-2600 
(bottom left) and PSR J1913— 0440 (bottom right). See caption for Fig. 2 and text for details. The choice of longitude zero for each 
pulsar is given in the text and is marked with a dashed line where an accurate determination can be made. 
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Figure 3. (u)-(x). Polarization profiles at 1.4 GHz for PSR J1921+2153 (top left), PSR J1932+1059 (top right), PSR J1935+1616 
(bottom left) and PSR J1941— 2602 (bottom right). See caption for Fig. 2 and text for details. The choice of longitude zero for each 
pulsar is given in the text and is marked with a dashed line where an accurate determination can be made. 
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little profile evolution with frequency and also shows a steep 
swing of polarization through the middle of the pulse. We 
have therefore located 0o at the symmetry centre of the pro- 
file. An RVM fit can be attempted for this pulsar; although 
the fit is not very good, the fitted <po is coincident with our 
chosen 4>o- 

PSR J0908-1739 (Figure 3f): This pulsar is clas- 
sified as a partial cone (likely the leading edge) by Lyne & 
Manchester (1988). There is virtually no frequency evolu- 
tion of the profile between 0.4 and 3.1 GHz lending support 
to this idea. In Figure 3f we have simply labelled the profile 
midpoint as longitude zero, although its location is likely to 
be at significantly later longitudes. We do not assign phio 
or PAo for this pulsar. 

PSR J0922+0638 (Figure 3g): This pulsar is classified 

as a partial cone (likely the trailing part) by Lync & Manch- 
ester (1988). Weisbcrg ct al. (1999) provide a comprehensive 
review of the many observations of this pulsar. It is likely 
that 4>o is located somewhat earlier than the centre of emis- 
sion, where the polarized emission is too low to determine 
PAo. 

PSR J0953-I-0755 (Figure 3h): The profile consists of a 
main and interpulse and has (polarized) emission over vir- 
tually the whole 360° of longitude. RYM fitting has been 
attempted by a number of authors as summarised in Everett 
& Weisberg (2001). The fit given by those authors shows a 
preference for a single-pole interpretation of the polaxization 
PA swing in this pulsar, unlike other authors who prefer an 
orthogonal rotator (e.g. Blaskiewicz, Cordcs & Wasserman 
1991). As Everett & Weisberg (2001) appear to have the 
best data, wc have used their fit to establish the location of 
4>o at a longitude which leads the peak of the main pulse by 
12°. 

PSR J1136-I-1551 (Figure 3i): The pulsar appears to 
have a classic double-cone profile and we locate 4'o at the 
midpoint of the profile. 

PSR J1239-I-2453 (Figure 3j): This pulsar is a classic 
five-component pulsar. It is likely that the core component 
is located close to the centre of the profile flanked by the 
inner and outer cones. This interpretation is backed up by 
the fact that at high frequencies the core component is less 
prominent and the swing of the handedness of circular polar- 
ization coincides with 4>o- The PA swing in this pulsar does 
not conform to the RVM so no value of 4>o can be obtained 
from an RVM fit. 

PSR J1430-6626 (Figure 3k): This pulsar has an ex- 
tremely complex PA curve at 1.4 GHz and appears to swing 
through 270° , not allowed in the RVM model. There is rather 
little profile evolution between 0.4 and 1.4 GHz but the 3.1 
GHz observations show a strong decline in the flux density 
of the leading component. At the peak of the trailing compo- 
nent, the sign of circular polaxization changes from negative 
to positive. It seems likely that 4>o is close to the pealc of this 
component, so we have used this location to estimate PAo. 
PSR J1453-6413 (Figure 31): The profile of this pul- 
sar evolves strongly with frequency. At high frequencies, the 
component leading the main peak dominates the profile and 
therefore seems likely to be a conal component. The promi- 
nent peak at 1.4 GHz is likely to be the core component and 
this component also shows a steep swing of PA. We therefore 
assign 0o to the peak of this component which also corre- 



sponds to a local infiexion point in the PA swing, although 
a formal fit to the RVM is not possible. 
PSR J1456-6843 (Figure 3m): Although there appear 
to be multiple components in the profile of this pulsar there 
is rather little evolution with frequency across a wide fre- 
quency range. The PA swing appears to be simpler at lower 
frequencies (Hamilton et al. 1977a), although this may be 
an artifact of low time resolution. We have placed (po where 
the sign of circular polarization changes; this is almost the 
centre of symmetry of the profile. 

PSR 1645-0317 (Figure 3n): The pulsar is classified 
as a core single at low frequencies as it shows a simple 
single-component profile (Hamilton et al. 1977a). At 5 GHz 
and above, the conal outriders become very prominent (von 
Hoensbroech & Xilouris 1997) and we are therefore confident 
that the peak of the profile marks 4>o- 

PSR J1709-1640 (Figure 3o): This pulsar is classified 

as a core single. Wc locate (fio at the peak of the profile, 
which is also the location at which the circular polarization 
changes sign. We were unable to estimate a value for the 
RM for this pulsar with sufficient accuracy to improve on 
the previously listed value. 

PSR J1709-4429 (Figure 3p): Unfortunately, polariza- 
tion observations of this pulsar are not particularly illumi- 
nating as it is difficult to tell whether <f>o is located near to 
the single, rather broad component or whether it is instead 
far out in the pulsar beam (as claimed by Manchester 1996 
and Crawford ct al. 2001). We therefore do no assign (f>o or 
PAo for this pulsar. 

PSR J1740+1311 (Figure 3q): We did not obtain high 

frequency observations of this pulsar and used the value 
of RM given in Weisberg et al. (2004). This pulsar is a 5- 
component pulsar with the core located close to the centre 
of the profile (see the discussion in Weisberg et al. 2004) 
and we assign 0o to this location. The PA swing also shows 
an inflection point at this same longitude, further strength- 
ening the case that this location marks the magnetic pole 
crossing. 

PSR J1844+1454 (Figure 3r): This pulsar only has a 
single component at low frequencies. At higher frequencies a 

conal outrider becomes prominent and dominates the profile 
above 1.4 GHz. We have located (po at the centre of the 
strong low frequency component. 

PSR J1900-2600 (Figure 3s): The profile is complex 
and shows multiple components. The circular polarization 

shows a characteristic swing with a change of handedness 
near the centre of symmetry of the profile. This fact, cou- 
pled with observations at lower frequencies (McCuUoch et al. 
1978) make it likely that (t>o is near the centre of the profile. 
PSR J 1913-0440 (Figure 3t): The profile consists of 
three components at 1.4 GHz. At lower frequencies the 
trailing component dominates the profile (Hamilton et al. 
1977a), whereas at higher frequencies the leading two com- 
ponents become more dominant. Given that there is also a 
steep swing of PA associated with the trailing component 
wc locate 0)0 near the peak of this component. 
PSR J1921-I-2153 (Figure 3u): Weisberg et al. (1999) 
describe observations of this pulsar made at various frequen- 
cies. The pulsar is likely a conal double with (f>() near the 
profile centre. The PA profile deviates strongly from that 
expected in the RVM model. 

PSR J1932+1059 (Figure 3v): This pulsar has emission 
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Table 2. PA„, PAq and their offset, |*|, for our sample of 25 
pulsars. The figure in brackets gives the error in the last digit(s). 



Jname 


Bname 


PA„ 


PAo 


1*1 


Fig 






(deg) 


(deg) 


(deg) 






Pulsars with I'l' <10° 


or 1*1 >80° 






J0630-2834 


B0628-28 


294(3) 


26(2) 


88(4) 


3c 


J0742-2822 


B0740-28 


278(5) 


-81.7(1) 


0(5) 


3d 


J0820-1350 


B0818-13 


159(6) 


65(2) 


86(6) 


3e 


J0835-4510 


B0833-45 


301.0(1) 


36.8(1) 


84.2(2) 


2 


J1239+2453 


B1237+25 


295.0(1) 


-66(1) 


1(1) 


3j 


J1430-6623 


B1426-66 


236(9) 


-28.5(7) 


85(9) 


3k 


J1453-6413 


B1449-64 


217(3) 


-56.9(4) 


86(3) 


31 


J1709-1640 


B1706-16 


192(16) 


15(2) 


3(16) 


3o 


J1740+1311 


B4737+13 


227(6) 


-46(4) 


87(7) 


3q 


J1844+1454 


B1842+14 


36(15) 


-52(2) 


88(15) 


3r 




Pulsars 


with 10° < 


1*1 < 80° 






J0525+1115 


B0523+11 


132(16) 


-65(4) 


17(16) 


3a 


J0953+0755 


B0950+08 


355.9(2) 


14.9(1) 


19.0(2) 


3h 


J1136+1551 


B1133+16 


348.6(1) 


-78(2) 


67(2) 


31 


J1456-6843 


B1451-68 


252.7(6) 


-31.6(6) 


76(1) 


3m 


J1645-0317 


B1642-03 


353(3) 


56(4) 


63(5) 


3n 


J1900-2600 


B1857-26 


202.8(7) 


-43(2) 


66(2) 


3s 


J1913-0440 


B1911-04 


166(11) 


-68(2) 


54(11) 


3t 


J1921+2153 


B1919+21 


34(12) 


-35(2) 


69(12) 


3u 


J1932+1059 


B1929+10 


65.2(2) 


-11.3(1) 


76.5(2) 


3v 


J1935+1616 


B1933+16 


176(1) 


10.1(7) 


14(1) 


3w 


Pulsars for which PAq 


and hence 


* cannot be determined 


J0601-0527 


B0559-05 


194(16) 






3b 


J0908-1739 


B0906-17 


167(2) 






3f 


J0922+0638 


B0949+06 


12.0(1) 






3g 


J1709-4429 


B 1706 -44 


160(10) 






3p 


J19 11-2(102 


L519:i7-2(i 


l:-!0(17) 









over a large longitude range and also has a low-amplitude 
interpulse; in the figure we show only the main pulse. For 
a more sensitive observation at this frequency see Weisberg 
ct al. (1999). Many authors have attempted an RVM fit for 
this pulsar, and these are discussed by Everett & Weisberg 
(2001). We use their fitted results to locate 4>o- 
PSR J1935+1616 (Figure 3w): This pulsar is classified 
as a core single at low frequencies. At high frequencies, conal 
outriders appear but the central part of the profile remains 
constant (see the discussion in Weisberg et al. 1999). The 
PA is rather complicated with a number of orthogonal jumps 
across the profile and an RVM fit is not possible. We there- 
fore take the location where the circular polarization changes 
hand to be (pQ. 

PSR J1941-2602 (Figure 3x): The profile of this pulsar 
also seems likely to be a partial cone. The PA variation is 

rather flat and there is virtually no profile evolution over a 
wide frequency range. In particular, the shoulder seen at the 
trailing edge of the profile remains constant with frequency. 
We cannot assign 4>o or PAo for this pulsar. 
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Figure 4. Histograms of observed (top) and simulated (bottom) 
values of I^I"]. The left hand panels ignore the possibility of or- 
thogonal mode emission so that PAo=PAr. The right hand panels 
include this possibility which then reduces the maximum offset to 
45° . Also included is measurement uncertainties in the data with 
typical value of 10°. Top panels: Observed values of in bins of 
10° for 20 pulsars. Bottom panels: Simulated histogram for PA^ 
parallel to PA„ (thick line), uncorrclated PA^ and PA^ (dashed 
line) and PA^ perpendicular to PA„ (dot-dash line). 

from Table 1. Column 5 shows the offset, \'^\, between PA„ 
and PAo where ^ is defined as 

* = PA„ - PAo; -90° < * < 90° (3) 

We force ^' to lie between —90° and 90° by rotating PAo by 
±180° as necessary (note that any polarized PA has by na- 
ture a 180° ambiguity). The error in \I' comes from a quadra- 
ture sum of the errors in PA„ and PAq. 

The table is subdivided into three sections. The top sec- 
tion lists those pulsars for which we find a good relationship 
between PAo and PA„ defined as where l^"] is less than 10° 
or greater than 80°. There arc 10 pulsars in this category. 
The second section lists the 10 pulsars for which there is no 
clear relation between PAo and PA^ i.e. where \^\ lies be- 
tween 10 and 80° . Section 3 lists the 5 pulsars where we can 
not make a judgement either because there is either no emis- 
sion at 00 or the PA swing is too complicated to determine 
PAo. 



6 DISCUSSION 

In Table 2 we summarise our results. Columns 3 and 4 show 
the position angles of the velocity vector, PA„, and of the 
linear polarisation, PAo, at 0o, the closest approach of the 
line-of-sight to the magnetic pole (see Equation 1), taken 



6.1 Statistics 

The top left panel of Figure 4 shows the observed distri- 
bution of I "If I as a histogram for the 20 pulsars for which 
we have measured If we assume that PA„ and PAo are 
unrelated then the expectation is that the offsets between 
them should be randomly distributed in the range 0° to 90° . 
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This is shown by the dashed line on the bottom left panel of 
Fig. 4. The probability of having 7 objects in the top bin, as 
observed, is very small, 2.6 x 10~*. A Kolmogorov-Smirnov 
(KS) test rules out a random distribution at the 94 per cent 
confidence level. Clearly our measurements show that PA^ 
and PAo must be related in some fcishion. Consider now 
some possible types of relationship between the observables 
PA„ and PAo; namely that they are either parallel or per- 
pendicular. To start, we ignore the confounding factor of the 
orthogonal emission so that PA()=PAr (sec Fig. 1). 

The case where PA^ is parallel to PA,, is a simple one - 
all the pulsars should fall in the first bin of the histogram. 
However, measurement errors will broaden this distribution 
as shown by the solid line on the bottom left panel of Fig. 4. 
This also does not represent the data; the KS test rules out 
that these distributions are the same at greater than the 
99.9 per cent confidence limit. 

The case where PAt,, is perpendicular to PAr is slightly 
more complex. In this case, the projection effects from the 
true 3-D vector to the 2-D sky vector will cause the (ob- 
served) ofTset between PA„ and PAo to spread away from 
90°. This is shown by the dash-dot line on the bottom left 
panel of Fig. 4 The result looks remarkably similar to the 
real data apart from the excess in the first two bins. The KS 
test here does not rule out the possibility that the distribu- 
tions are the same. 

However, the statistics do not take into account the pos- 
sible presence of orthogonal mode omission which compli- 
cates the underlying relationship between PA^ and PAo, as 
shown in Fig. 1. Either PA^ is parallel to PAo if the pulsar 
is emitting in the 'normal' mode or PAr is perpendicular to 
PAo if the pulsar emission is in the 'orthogonal mode'. Si- 
multaneously allowing for either possibility serves to reduce 
the maximum offset between PAo and PA„ to only 45°. The 
right hand panel of Fig. 4 shows the data and the simulations 
for such a scenario. The observed histogram appears some- 
what broader than the sinmlation but this may be caused 
by the possible sources of contamination discussed below. 
In any case, the results show that we can confidently rule 
out that PA,, and PA„ are unrelated (dashed line) at the 
98 per cent confidence limit. The KS test also rules out the 
case where PAr is perpendicular to PA„ (dash-dot line) at 
greater than the 99 per cent confidence limit but does not 
rule out that the case where PA,, is parallel to PA^ (solid 
line) is the same as the observed distribution. 

6.2 Information from outside the radio-band 

In the previous section we have demonstrated that our ra- 
dio observations provide clear evidence for a relationship be- 
tween the spin and velocity vectors of pulsars. The remain- 
ing uncertainty as to whether parallel and orthogonal vec- 
tors occur or whether the orthogonal PA values are caused 
by orthogonal emission modes, can be decided by using in- 
formation obtained outside the radio window. We consider 
these additional pieces of information in turn. 

6.2.1 The Vela pulsar 

In Section 1 we noted that, in the Vela pulsar, the axis of 
the X-ray torus (and hence the pulsar's rotation axis) had a 



position angle on the sky of 130° and 310° and that this com- 
pared well with the proper motion position angle of 301°. 
In order to compare this with the polarization position an- 
gle, Helfand et al. (2001) used the results of Hamilton et 
al. (1977) who measured the polarization PA at the peak of 
the total intensity and obtained an angle of 64° (or 244°). 
This angle is about 60° different from the position angle 
of the torus and proper motion. Helfand et al. (2001) then 
make the bold assertion that as the angle is "almost" 90°, 
the emission from Vela must be predominantly in a mode 
polarized orthogonally to the magnetic field lines. Radhakr- 
ishnan & Deshpande (2001) discussed the implications of 
the Vela pulsar's orthogonal mode emission in the context 
of the physical model of Luo & Melrose (1995). 

There are two issues here. First, 60° is not that close to 
90°. Secondly, the peak emission of the profile almost cer- 
tainly does not coincide with <po. As described above, this is 
best measured at the point of greatest swing of the PA curve 
which, in the Vela pulsar, occurs significantly later than the 
pulse peak. In our data, the PA at this location appears to 
be almost exactly 90° offset from PA„. This provides the 
observational evidence that the hypothesis by Helfand et al. 
(2001) and Radhakrishnan & Deshpande (2001) is correct 
and it is likely that the emission in Vela is perpendicular to 
the magnetic field direction and that PA„ and the rotation 
axis are in fact aligned. 

6.2.2 The Crab pulsar 

Two pieces of evidence from the Crab pulsar also point 
in this direction. The Crab has PA„ of 292° (Caravoo & 
Mignani 1999), and interpretation of the X-ray torus indi- 
cates that PA,, in the Crab is 124° or 304° (Ng & Romani 
2004) meaning that PA^ is parallel to PA„ within the mear- 
surement uncertainties. 

Secondly, while the radio polarization of this pulsar is 
very complex, polarization measurements in the optical and 
UV show a much simpler picture. The main and inter pulses 
show PA variations which are consistent with a RVM and 
both have a PAo near 124° (Smith et al. 1988; Graham-Smith 
et al. 1996), very similar to the PA of the X-ray torus. In 
the radio, at low frequencies the PA of the main and in- 
terpulses is 120°, whereas at higher frequencies the PA of 
the interpulse is 30° (Moffett & Hankins 1999). Although it 
is unclear where the main and interpulse emission actually 
arises (there are competing models for polar cap and light 
cylinder emission) , it is interesting that the polarization PAs 
match both PA„ and PA^. 



6.2.3 Other pulsars 

PSR B0656-hl4 has PA„ of 93° (Brisken et al. 2003b), and 

recent optical polarization measurements by Kern et al. 
(2003) match well with the RVM fit of Everett & Weisberg 
(2001) allowing a determination of PAq of 98°. As we do 
not have to worry about orthogonal mode emission at opti- 
cal wavelengths, then PA, is parallel to PA„. For the young 
pulsar PSR B0540-69 in the Large Magellanic Cloud, Scr- 
afimovich et al. (2004) have shown (weak) evidence that the 
velocity vector and the jet (and hence rotation axis) align. 
Finally, other pulsars in the sample of Ng & Romani (2004), 
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PSRs J0538+2817, B1951+32 and B1706-44 also show evi- 
dence of tori which also indicate that PA,, is parallel to PA„. 

6.3 Possible contaminants 

We are aware of a number of effects that may influence 
the observed PA„ and PAo values, and which may serve 
to broaden the observed distribution of >]/. Apart from mea- 
surement uncertainties and orthogonal mode emission al- 
ready discussed, we can identify four possible sources of 
contamination. These are (i) propagation of the radiation 
through the pulsar magnetosphere, (ii) the effects of differ- 
ential galactic rotation, (iii) the effects of the galactic po- 
tential on the pulsar's peculiar motion and (iv) the motion 
of the pre-supernova progenitor star. We will deal with each 
of these in turn. 

6.3.1 Magnetospheric propagation 

Blaskiewicz et al. (1991) considered the effects of retarda- 
tion and aberration of radiation travelling through the pul- 
sar magnetosphere. In brief their argument implies that the 
position angle sweep lags the profile by an amount which 
depends on the altitude at which the emission is radiated. A 
consequence of this is that if 4>o is determined from the RVM 
fit, it should lag the 0o determined from the profile char- 
acteristics. Furthermore, under the paradigm that different 
frequencies are emitted from different heights in the magne- 
tosphere the effect should be frequency dependent. We have 
estimated the size of the effect by computing an emission 
height at our observing frequency using the prescription in 
Kijak & Gil (2003) and then deriving the lag between the 
PA curve and the intensity profile. We find that in the ma- 
jority of cases this effect is loss than 3°, but for the Vela 
pulsar and PSR J0742-2822 the lag is 13.5 and 8.5° respec- 
tively. As we have demonstrated, other pieces of evidence 
suggest that for Vela our determination of is correct and, 
in any case, for Vela we have directly measured (po from the 
RVM and hence it does not depend on the aberration value. 
For PSR J0742-2822, Karastergiou & Johnston (2005) have 
shown that the lag must be significantly smaller than that 
given in Kijak & Gil (2003). As the derivation of the emis- 
sion heights of pulsars is highly model dependent we have 
not included these values in our error budget. 

6.3.2 Differential galactic rotation 

The proper motion vector of the pulsar at birth is what we 
really want to compare with PA,-, since the latter quantity 
will remain fixed in the absence of further torques. Prior to 
birth, the pulsar's progenitor star was orbiting the Galac- 
tic Centre with a velocity which is location dependent. The 
difference between this velocity and the orbital velocity at 
the Sun's location is the differential galactic rotation com- 
ponent. The measurement of the proper motion is made 
with respect to the solar system barycentre and therefore 
includes the effects of differential galactic rotation on the 
measured value. Since birth, the pulsar has moved in the 
Galactic potential which has therefore changed its velocity 
and proper motion vector. In order to determine the rotar 
tion of the proper motion vector over the pulsar's lifetime, it 



is necessary to integrate equations of motion in the Galac- 
tic potential, using as boundary conditions the pulsar's age, 
distance, and three-dimensional velocity. The pulsar's char- 
acteristic age, Tc, listed in Table 1 gives a useful indication 
of the true age but may be in error by a factor of 2 and pos- 
sibly significantly more for young pulsars [e.g. in Vela (Lyne 
et al. 1996) and PSR J0538-I-2817 (Kramer et al. 2003)]. The 
pulsar distance can be determined from the parallax (where 
known) or by converting the dispersion measure using the 
model of Cordes & Lazio (2005) which has an uncertainty for 
any given object of up to 50 per cent. The observed proper 
motion gives two components of the velocity, but the radial 
component is completely unknown although one can justify 
that it must lie in the range ±1000 kms^^. A possible handle 
oil the radial velocity can be obtained in the following way. 
If we assume that the pulsar was born in the Galactic plane 
and we know its current location and age today, we can ad- 
just the radial velocity so that extrapolating backwards in 
time forces a birth place in the Galactic plane. Such a tech- 
nique has boon used relatively successfully in the case of the 
relativistic binary pulsars (Wex, Kalogera & Kramer 2000; 
Willems, Kalogera & Henninger 2004). However, given the 
large uncertainties in the age and distance to the pulsars, 
we do not feel this approach is appropriate here. 

Estimating the differential galactic rotation effect on 
PAy at birth is therefore difficult as we do not know the 
pulsar's birth location. We can make some general state- 
ments, however, using the potential model of Kuijken & 
Gilmore (1989) which yields a rotation speed of the Sun of 
221.4 kms~^. For young pulsars, born within ~1.5 kpc of the 
Sun, the differential galactic rotation contributes no more 
than 10-20 kms"'^ to the overall velocity. We have simulated 
this effect on our sample and find that it will contribute no 
more than 5° to the error budget except for very slow mov- 
ing pulsars. For the most distant pulsar in our sample, PSR 
J1935-I-1616 (for which ^'=14°), at a distance in excess of 
~5 kpc, the effect could contribute as much as 20° to the er- 
ror bar. It is therefore important to ensure that the sample 
contains young, nearby pulsars to minimise this effect. 

6.3.3 Acceleration in the galactic gravitational potential 

The effect of the galactic potential on the pulsar's pecu- 
liar motion depends critically on the unknown radial veloc- 
ity. Again, we can make only general statements about the 
magnitude of this effect. For pulsars born in the vicinity of 
the Sun, and moving with a velocity of ~300 kms^^, the 
time-scale for the proper motion vector to have significantly 
changed is 15 Myr. We note that 4 of the 5 pulsars with 
characteristic ages of 3 Myr or less have either |$| <10° 
or l^*! >80°. Conversely, of the 6 pulsars with ages greater 
than 10 Myr, only one (PSR J1239+2453) has close match 
between PAo and PA„. For the 9 remaining pulsars, with 
intermediate ages, 5 show a correlation between PAo and 
PA„. 



6.3.4 Motion of the progenitor star 

Finally, an important effect can occur when the pre- 
supernova star already has a peculiar velocity. Many high- 
mass stars are in binary pairs often with short orbital peri- 
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ods. The pre-supcrnova star can therefore aheady have a ve- 
locity of ^100 kms~^ or more. If the imparted kick velocity 
is also only of this order, then the resultant velocity vector 
will not align with the rotation axis. This scenario weis al- 
ready discussed by Doshpando ct al. (1999) who considered 
it a major reason for their apparent lack of correlation. 

6.4 Summary 

In summary, we have shown that the case for an alignment 
of the velocity and spin vectors is very strong. The case 
is strengthened by the fact that the youngest pulsars in the 
sample, which suffer from the smallest unknown errors, show 
the best alignment. In addition, we know that in the case 
of Vela, the observed emission is in an orthogonal mode. 
Wc think it likely therefore that PAi, and PAr are parallel 
in all cases but that many (though not all) pulsars emit 
predominantly in the orthogonal mode. This fits with the 
values returned from the KS tests on the histrograms in 
Fig 4 as discussed in section 6.1. 



7 IMPLICATIONS 

We see no compelling reason to revive the old rocket model 
of Tademaru & Harrison (1975) even though that model 
predicts that PA„ should be parallel to PAr. The prob- 
lems with that model still remain; it requires pulsar spins of 
1 ms at birth, it has difHculty explaining the binary neutron 
star systems, and polarization data do not favour off-centred 
dipoles. 

In the context of the Spruit & Phinney (1998) model 
the parallel and orthogonal cases arise in very different ways. 
In the perpendicular case, the number of impulses must be 
small and of short duration. In the parallel case, the impulses 
must be of long duration and can be significant in number. 
In principle, the two cases lead to different observables. In 
the perpendicular case the (initial) spin period and veloc- 
ity should be correlated. In the parallel case, short period 
pulsars should have a lower mean velocity. In practise, how- 
ever, distinguishing these cases observationally is difficult, 
although we note that there is some evidence that short pe- 
riod pulsars are slow moving. Theoretical arguments based 
on neutrino scattering and parity violation tend to favour 
long duration kicks but this might only be applicable in 
magnetic fields greater than 10^^ G (Arras & Lai 1999), 
values which seem unlikely for the radio pulsars discussed 
here. More recently, Schmitt, Shovkovy & Wang (2005) have 
speculated that a neutrino propulsion mechanism associated 
with the presence of a colour super-conductor in the trans- 
verse A phase, produces both high kick velocities and aligned 
spin and velocity vectors. It is clear that there are theoret- 
ical possibilities still to be explored (see also Lai, Chernoff 
& Cordos 2001). 

Both Iben & Tutokov (1996) and Spruit & Phinney 
(1998) argue that pre-supernova stars do not possess suf- 
ficient angular momentum to explain the fast spins of new- 
born neutron stars. This leads the latter authors to propose 
their model where one or more off-centre natal kicks is the 
principal source of young pulsars' short periods. We wish to 
point out an important general consequence of such models 
which seems to have escaped notice, namely that off-centre 



kick(s) at a radius r that significantly change the magnitude 
of the pulsar's spin angular momentum will probably also 
change the angular momentum direction. This is clear from 
the vector nature of the equation governing changes in angu- 
lar momentum, dL/dt — rx F, whore F is the force applied. 
Unless off-centre kick(s) are somehow applied in a fashion 
that preserves the original rotational symmetry of the star, 
the angular momentum vector will tip significantly. 

Observations of the current parameters of the dou- 
ble neutron star systems allow a determination of the 
pre-supernova evolution of the system (Wex, Kalogera 
& Kramer 2000; Willems, Kalogera & Henninger 2004; 
Thorsctt, Dewey & Stairs 2005). In turn, the kick magni- 
tude and direction imparted to the (second-born) neutron 
star can be constrained. In PSRs B1913+16, B1534-I-12 and 
J0737— 3039 the kick direction is significantly different from 
the spin axis of the pre-supernova star and the orbital an- 
gular momentum, with the inference that aligned kicks are 
therefore ruled out. However if, as shown above, the kick 
imparts both linear and angular momentum to the neutron 
star, then the spin axis of the newly-born neutron star is 
not related to either that of its progenitor, or the orbital 
angular momentum. Consequently, the derived constraints 
on the kick direction can shed no light on the alignment of 
spin and velocity vectors. 

There is an observational test which can be made. If 
Spruit & Phinney (1998) are correct, then the misalign- 
ment between the spin and orbital angular momenta in the 
PSR J0737-3039 system should be different for the A and B 
pulsars whereas if the kick only imparts velocity and leaves 
the pre-supernova spin axis unchanged then the misalign- 
ment angles should be identical. These measurements should 
be possible in the near future. Note that even if the misalign- 
ment angles are the same, the spin axes of the two pulsars 
will not be parallel at the present epoch since their preces- 
sion rates are different. 



8 CONCLUSIONS 

We have observed 25 pulsars at 1369 MHz, 21 of which were 
also observed at 3100 MHz. Careful polarization calibra- 
tion, including accurate RM determination, has enabled us 
to compute the PA of the polarized radiation at the pulsar 
to within an accuracy of 2-3°. We have used all available 
information to deduce <j)o, the longitude of the closest ap- 
proach to the magnetic pole of the star. We then compared 
the PA of the radiation at this longitude, which is related to 
the position angle of the axis of rotation, with the velocity 
vector derived from proper motion measurements. 

The combination of precise polarization measurements, 
accurate RMs and improved velocity information makes this 
study much more sensitive to any underlying correlation be- 
tween the spin and velocity vectors. In contrast to previous 
studies we restricted our sample to young, nearby pulsars 
where the effects of contaminants discussed earlier have lit- 
tle or no impact. 

We find a clear relationship between the velocity vec- 
tor and the rotation axis. The presence of orthogonal modes 
in pulsar emission makes it difficult to determine whether 
the velocity vector is parallel to the rotation axis or per- 
pendicular to it. However, additional information from the 
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optical and X-ray bands for PSR B0656+14, the Crab and 
Vela pulsars allows us to break this ambiguity: the velocity 
vector is parallel to the rotation axis and many pulsars emit 
linear polarisation with PA predominantly perpendicular to 
the magnetic field lines. 

We point out that if the kicks impart both linear and 
angular momentum to the new-born neutron star, the in- 
formation about the stellar axis prior to the supernova ex- 
plosion is lost. Observations of double neutron star systems 
cannot therefore rule out aligned kicks, contrary to claims 
in the literature. 
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